Adaptive tuning of Majorana fermions in a quantum dot chain. 
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We suggest a way to overcome the obstacles that disorder and high density of states pose to the 
creation of unpaired Majorana fermions in one-dimensional systems. This is achieved by splitting 
the system into a chain of quantum dots, which are then tuned to the conditions under which the 
chain can be viewed as an effective Kitaev model, so that it is in a robust topological phase with 
well-localized Majorana states in the outermost dots. The tuning algorithm that we develop involves 
controlling the gate voltages and the superconducting phases. Resonant Andreev spectroscopy allows 
us to make the tuning adaptive, so that each pair of dots may be tuned independently of the other. 
The calculated quantized zero bias conductance serves then as a natural proof of the topological 
nature of the tuned phase. 

PACS numbers: 74.45. +c, 74.78. Na, 73.63.Kv, 03.65.Vf 



I. INTRODUCTION 

Majorana fermions are the simplest q uas iparticles pre- 
dicted to have non-Abelian statistics!^ These topo- 
logically protected states can be realized in condensed 
matter systems, by making use of a combination of 
strong spin-orbit coupling, superconductivity and bro- 
ken time-reversal symmetry! 3 ^ Recently, a series of ex- 
periments have reported the possible observation of Ma- 
jorana fermions in semiconducting nanowiresf^l attract- 
ing much attention in the condensed matter community. 

Associating the observed experimental signatures ex- 
clusively with these non-Abelian quasiparticles, however, 
is not trivial. The most straightforw ard s ignature, the 
zero bias peak in Andreev conductanc e ! 10 ^ 1 ^ is not unique 
to Majorana fermions, but can appear as a result of var- 
ious physical mechanisms E2HHI such as the Kondo effect 
or weak anti-localization. It has also been pointed out 
that disorder has a detrimental effect on the robustness 
of the topological phase, since in the absence of time- 
reversal symmetry it may close the induced supercon- 
ducting gapP^This requires experiments performed with 
very clean systems. Additionally, the presence of multi- 
ple transmitting modes reduces the amount of control one 
has over such systemsp^U an( j ^he contribution of extra 
modes to conductance hinders the observation of Majo- 
rana fermions P3l Thus, nanowire experiments need setups 
in which only few modes contribute to conductance. 

In this work we approach the problem of realizing 
systems in a non-trivial topological phase from a dif- 
ferent angle. We wish to emulate the Kitaev chain 
modeP^ which is the simplest model exhibiting unpaired 
Majorana bound states. The proposed system consists 
of a chain of quantum dots (QDs) defined in a two- 
dimensional electron gas (2DEG) with spin orbit cou- 
pling, in proximity to superconductors and subjected to 
an external magnetic field. Our geometry enables us to 
control the parameters of the system to a great extent 
by varying gate potentials and superconducting phases. 




FIG. 1: Examples of systems allowing implementation of a 
Kitaev chain. Panel (a): a chain of quantum dots in a 2DEG. 
The QDs are connected to each other, and to superconductors 
(labeled SC), by means of quantum point contacts. The first 
and the last dots are also coupled to external leads. The 
normal state conductance of quantum point contacts (QPCs) 
between adjacent dots or between the end dots and the leads 
is G||, and of the QPCs linking a dot to a superconductor is 
G± . The confinement energy inside each QD can be controlled 
by varying the potential Vgate- Panel (b): Realization of the 
same setup using a nanowire, with the difference that each 
dot is coupled to two superconductors in order to control the 
strength of the superconducting proximity effect without the 
use of QPCs. 



We will show how to fine tune the system to the so- 
called "sweet spot" in parameter space, where the Ma- 
jorana fermions are well-localized at the ends of the sys- 
tem, making the topological phase maximally robust. A 
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sketch of our proposed setup is presented in Fig. . 

The setup we propose and the tuning algorithm is not 
restricted solely to systems created in a two-dimensional 
electron gas. The essential components are the ability 
to form a chain of quantum dots and tune each dot sep- 
arately. In semiconducting nanowires the dots can be 
formed from wire segments separated by gate-controlled 
tunnel barriers, and all the tuning can be done by gates, 
except for the coupling to a superconductor. This cou- 
pling, in turn, can be controlled by coupling two super- 
conductors to each dot and applying a phase difference 
to these superconductors. The layout of a nanowire im- 
plementation of our proposal is shown in Fig. [TJs). 

This geometry has the advantage of eliminating many 
of the problems mentioned above. By using single level 
quantum dots, and also quantum point contacts (QPC) 
in the tunneling regime, we solve issues related to mul- 
tiple transmitting modes. Additional problems, such as 
accidental closings of the induced superconducting gap 
due to disorder, are solved because our setup allows us 
to tune the system to a point where the topological phase 
is most robust, as we will show. 

We present a step-by-step tuning procedure which fol- 
lows the behavior of the system in parallel to that ex- 
pected for the Kitaev chain. As feedback required to 
control every step we use the resonant Andreev conduc- 
tance, which allows to track the evolution of the system's 
energy levels. We expect that the step-by-step structure 
of the tuning algorithm should eliminate the large num- 
ber of non-Majorana explanations of the zero bias peaks. 

A related layout together with the idea of simulating 
a Kitaev chain was proposed recently by J. D. Sau and 
S. Das SarmaP^ Although similar in nature, the geom- 
etry which we consider has several advantages. First of 
all, coupling the superconductors to the quantum dots in 
parallel, allows us to not rely on crossed Andreev reflec- 
tion. More importantly, being able to control inter-dot 
coupling separately from all the other properties allows 
to address each dot or each segment of the chain elec- 
trically. This can be achieved by opening all the QPCs 
except for the ones that contact the desired dots. 

This setup can also be extended to more complicated 
geometries which include T-junctions of such chains. 
Benefiting from the high tunability of the system and 
the localization of the Majorana f ermio ns, it might then 
be possible to implement braiding p 7 * 28 ! and demonstrate 
their non-Abelian nature. 

The rest of this work is organized as follows. In sec- 
tion |TT] we briefly review a generalized model of Kitaev 
chain, and identify the "sweet spot" in parameter space 
in which the Majorana fermions are the most localized. 
The system of coupled quantum dots is described in sec- 
tion [TTTJ For the purpose of making apparent the resem- 
blance of the system to the Kitaev chain, we present a 
simple model which treats each dot as having a single 
spinful level. We then come up with a detailed tuning 
procedure describing how one can control the parameters 
of the simple model, in order to bring it to the desired 



point in parameter space. In section |IV| our tuning pre- 
scription is applied to the suggested system of a chain of 
QDs defined in a 2DEG, and it is shown using numeri- 
cal simulations that at the end of the process the system 
is indeed in a robust topological phase. We conclude in 
section [V] 



II. GENERALIZED KITAEV CHAIN 

In order to realize unpaired Majorana bound states, 
we start from the Kitaev chairP^ generalized to the case 
where the on-site energies as well as the hopping terms 
are not uniform and can vary from site to site. The gen- 
eralized Kitaev chain Hamiltonian is defined as 

L-l 



n=l 



h.c. ) + e n a} n a n 



(2.1) 



where a n are fermion annihilation operators, e„ are the 
on-site energies of these fermions, t n cxp(i9 n ) are the 
hopping terms, and A„ exp(i(f> n ) are the p-wave pairing 
terms. 

The chain supports two Majorana bound states local- 
ized entirely on the first and the last sites, when (i): 
£„ = 0, (ii): A„ = t n , and (hi) (f> n+1 - <t> n - 6 n+1 - Q n = 0. 
The larger values of t n lead to a larger excitation gap. 
The condition (iii) is equivalent, up to a gauge transfor- 
mation, to the case where the hopping terms are all real, 
and the phases of the p-wave terms are uniform. The en- 
ergy gap separating the Majorana modes from the first 
excited state then equals 



-Egap = 2min{t n } r 



(2.2) 



The above conditions (i)-(iii), constitute the "sweet 
spot" in parameter space to which we would like to tune 
our system. Since all of these conditions are local and 
only involve one or two sites, our tuning procedure in- 
cludes isolating different parts of the system and mon- 
itoring their energy levels. For that future purpose we 
will use the expression for excitation energies of a chain 
of only two sites with e\ = £2 = 0: 



E 12 = ±{t 1 ±A 1 ). 



(2.3) 



III. 



SYSTEM DESCRIPTION AND THE 
TUNING ALGORITHM 



The most straightforward way to emulate the Kitaev 
chain is to create an array of spinful quantum dots, and 
apply a sufficiently strong Zeeman field such that only 
one spin state stays close to the Fermi level. Then the 
operators of these spin states span the basis of the Hilbert 
space of the Kitaev chain. If we require normal hopping 
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between the dots and do not utilize crossed Andreev re- 
flection, then in order to have both t n and A„ nonzero 
we need to break the particle number conservation and 
spin conservation. The former is achieved by coupling 
each dot to a superconductor, the latt er ca n be achieved 
by spatially varying Zeeman coupling J22I22I or more con- 
ventionally by using a material with a sufficiently strong 
spin-orbit coupling. Examples of implementation of such 
a chain of quantum dots in a two dimensional electron gas 
and in semiconducting nanowires are shown in Fig. [I] 

We neglect all the levels in the dots except for the one 
closest to the Fermi level, which is justified if the level 
spacing in the dot is larger than all the other Hamilto- 
nian terms. We neglect the Coulomb blockade, since we 
assume that the dot is strongly coupled to the supercon- 
ducting lead. The general form of the BdG Hamiltonian 
describing such a chain of spinful single-level dots is then 
given by: 

H s = ^ (Mn°o + V z a z ) c\ s c n ^ 

+ ^(A ind)n e««^ct )S ct iS ,+h.c.) (3 ' 1} 
+ {w n e lXn ' T c] l s c n+ i. s i + h.c.) , 

where c* n s and c rhS are creation and annihilation oper- 
ators of a fermion with spin s in the n-th dot, and Oi 
are Pauli matrices in spin space. The physical quanti- 
ties entering this Hamiltonian are the chemical poten- 
tial n n , the Zeeman energy V z , the proximity-induced 
pairing Ai n d, n exp(i$ n ), and the inter-dot hopping w n . 
The vector A„ characterizes the amount of spin rotation 
happening during a hopping between the two neighbor- 
ing dots (the spin rotates by a 2|A| angle). This term 
may be generated either by a spin-orbit coupling, or by 
a position-dependent spin rotation, required to make the 
Zeeman field point in the local z-directionPsH^H The in- 
duced pairing in each dot A; nc j in exp («<&„) is not to be 
confused with the p-wave pairing term A„exp(i0 n ) ap- 
pearing in the Kitaev chain Hamiltonian (2.1). 



In order for the dot chain to mimic the behavior of the 
Kitaev chain in the sweet spot, each dot should have a 
single fermion level with zero energy, so that e n = 0. Di- 
agonalizing a single dot Hamiltonian yields the condition 
for this to happen: 



ind. 



(3.2) 



When this condition is fulfilled, each dot has two 
fermionic excitations 

'Vz - A*n <4 t - e"'*" \/V z + n n Cnlj (3.3) 

bn = —7= [s/Vz ~ + e~'*" y/V z + fin C„f) . (3.4) 

The energy of a n is zero, the energy of b n is 2V Z . If the 
hopping is much smaller than the energy of the excited 




state, w n <C V z , we may project the Hamiltonian (3.1) 
onto the Hilbert space spanned by a n . The resulting 
projected Hamiltonian is identical to the Kitaev chain 
Hamiltonian of Eq. ( 2.1 ), with the following effective pa- 
rameters: 

e„ = 0, (3.5a) 
t n e l9n = w n (cos \ n + i sin X n cos p n ) x 

sin (a n+1 + a n ) cos(<5$„/2) (3.5b) 

+i cos (a„ + i - a n ) sin(<5$„/2) , 
A n e l< ^™ = iw n sin A„ sin p n e'^" x 

cos(a n+ i +a„)cos(5$„/2) (3.5c) 

+i sin (a„ + i - a n ) sin (<5$ n /2) 



where 



fi n = V z sin(2a„), A ilad ,n = V z cos(2a„), 



(3.6) 



K = A„ (sinp„cos^„, sinp„sin^„, cosp„) , (3.7) 

and 6$ n = $„ - $ n+ i. 

It is possible to extract most of the parameters of the 
dot Hamiltonian from level spectroscopy and then tune 
the effective Kitaev chain Hamiltonian to the sweet spot. 
The tuning, however, becomes much simpler if two out 
of three of the dot linear dimensions are much smaller 
than the spin-orbit coupling length. Then the direction of 
spin-orbit coupling does not depend on the dot number, 
and as long as the magnetic field is perpendicular to the 
spin-orbit field, the phase of the prefactors in Eqs. (3.5) 
becomes position-independent. Additionally if the dot 
size is not significantly larger than the spin-orbit length, 
the signs of these prefactors are constant. This ensures 
that if <5$ n = 0, the phase matching condition of the 
Kitaev chain is fulfilled. Since 8$> n = leads to both t n 
and A„ having a minimum or maximum as a function 
of (5$„, this point is straightforward to find. The only 
remaining condition, t n = A n at <5$ = 0, requires that 

a n + a n+l = Ki- 

The above calculation leads to the following tuning 
algorithm: 

1. Open all the QPCs, except for two contacting a 
single dot. By measuring conductance while tuning 
the gate voltage of a nearby gate, ensure that there 
is a resonant level at zero bias. After repeating for 
each dot the condition e„ = is fulfilled. 

2. Open all the QPCs except the ones near a pair of 
neighboring dots. Keeping the gate voltages tuned 
such that s n = 0, vary the phase difference between 
the neighboring superconductors until the lowest 
resonant level is at its minimum as a function of 
phase difference, and the next excited level at a 
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maximum. This ensures that the phase tuning con- 
dition 4> n +\ — 4>n — #n+i — On = is fulfilled. Repeat 
for every pair of neighboring dots. 

3. Start from one end of the chain, and isolate pairs 
of dots like in the previous step. In the pair of n-th 
and n + 1-st dots tune simultaneously the coupling 
of the n + 1-st dot to the superconductor and the 
chemical potential in this dot, such that e n +\ stays 
equal to 0. Find the values of these parameters such 
that a level at zero appears in two dots when they 
are coupled. After that proceed to the following 
pair. 

Having performed the above procedures, the coupling 
between all of the dots in the chain is resumed, at which 
point we expect the system to be in a robust topological 
phase, with two Majorana fermions located on the first 
and last dots. In practice one can also resume the cou- 
pling gradually by, for instance, isolating triplets of adja- 
cent dots, making sure they contain a zero-energy state, 
and making fine-tuning corrections if necessary, and so 
on. 



IV. TESTING THE TUNING PROCEDURE BY 
NUMERICAL SIMULATIONS 

We now test the tuning procedure by applying it to a 
numerical simulation of a chain of three QDs in a 2DEG. 
The two-dimensional BdG Hamiltonian describing the 
entire system of the QD chain reads: 

(p 2 \ a 

— + V(x, y)\T z + jipxVy - T z a y p x ) 

+ A ind (cos($)r a + sin($)r x ) a y + V z t z <t z . 

(4.1) 

Here, o~i and T; are Pauli matrices acting on the spin and 
particle-hole degrees of freedom respectively. The term 
V{x 1 y) describes both potential fluctuations due to dis- 
order, and the confinement potential introduced by the 
gates. The second term represents Rashba spin-orbit cou- 
pling, Ai n( j(a;, y) -exp ($(#, y)) is the s-wave superconduc- 
tivity induced by the coupled superconductors, and V z is 
the Zeeman splitting due to the magnetic field. Full de- 
scription of the tight-binding equations used in the sim- 
ulation is presented in Appendix |A") 

The chemical potential of the dot levels n n is tuned by 
changing the potential V(x, y). For simplicity we used 
a constant potential V n added to the disorder potential, 
such that V(x,y) = V n + Vo(x, y) in each dot. Varying 
the magnitude of Aj nc i.n is done by changing conductance 
G± of the quantum point contacts, which control the 
coupling between the dots and the superconductors. Fi- 
nally, varying the superconducting phase $(x, y) directly 
controls the parameter $„ of the dot to which the su- 
perconductor is coupled, although they need not be the 
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FIG. 2: Andreev conductance measured from the left lead as a 
function of bias voltage and QD potential (measured relative 
to quarter filling) for the second dot. Changing the chemical 
potential allows to tune quasi-bound states to zero energy 
(white circle). 



The tuning algorithm required monitoring the energy 
levels of different parts of the system. This can be 
achieved by measuring the resonant Andreev conduc- 
tance fro m one of the leads. The Andreev conductance 
is given b> |33 1 34 1 

G/G = N- Tr(r ee r|j + Tr^,^), (4.2) 

where Go = e 2 /h, N is the number of modes in a given 
lead, and r ee and r^ e are normal and Andreev reflection 
matrices. Accessing parts of the chain (such as a single 
dot or a pair of dots) can be done by opening all inter- 
dot QPCs, and closing all the ones between dots and 
superconductors, except for part of the system that is of 
interest. 

We begin by finding such widths of QPCs that Gii ~ 
0.02 and G± ~ 4Go- This ensures that conductance be- 
tween adjacent dots, is in tunneling regime and that the 
dots are strongly coupled to the superconductors such 
that the effect of Coulomb blockade is reduced.^ The 
detailed properties of QPCs are descrbed in App. [A|and 
their conductance is shown in Fig. [8] 

First step: tuning chemical potential. We sequentially 
isolate each dot, and change the dot potential V n . The 
Andreev conductance as a function of V n and bias voltage 
for the second dot is shown in Fig. [2j We tune V n to the 
value where a conductance resonance exists at zero bias. 
This is repeated for each of the dots and ensures that 
fi n = 0. 

Second step: tuning the superconducting phases. We 
now set the phases of the induced pairing potentials $„ 
to constant. As explained in the previous section, this oc- 
curs when A n and t n experience their maximal and min- 
imal values. According to Eq. ( |2.3| ) this happens when 
the separation between the energy levels of the pair of 
dots subsection is maximal. Fig. [3] shows the evolution 
of these levels as a function of the phase difference be- 
tween the two superconductors. The condition 5&i = 
is then satisfied at the point where their separation is 
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SC phase difference 



FIG. 3: Conductance as a function of bias voltage and super- 
conducting phase difference for a two-dot system. The two 
lowest energy levels are given by Eq. ( |2.3[ ) of a two site Ki- 
taev chain, as indicated. At the point where their separation 
is maximal (SC phase difference in the plot), the phase dif- 
ference <5$ n of the induced superconducting gaps vanishes. 



maximal. 

Third step: tuning the couplings. Finally we tune t n = 
A„. This is achieved by varying G_l, while tracking the 
Andreev conductance peak corresponding to the t n — A„ 
eigenvalue of the Kitaev chain we are emulating. After 
every change of Gj_ we readjust V n in order to make sure 
that the condition e n = (or equivalently V* z 2 = (j% + A 2 ) 
is maintained. This is necessary because not just A„, but 
also [i n depend on G±. Therefore, successive changes of 
G± and V n are performed until the smallest bias peak is 
located at zero bias. The tuning steps of the first two 
dots are shown in Fig. [4] We repeat steps 2 and 3 for 
each pair of dots in the system. 

Finally, having full all three conditions required for 
a robust topologically non-trivial phase, we probe the 
presence of localized Majorana bound state in the full 
three-dot system by measuring Andreev conductance (see 
Fig. [5]). In this specific case, the height of the zero bias 
peak is approximately 1.85 Go, signaling that the end 
states are well but not completely decoupled. Increasing 
the transparency of the QPC connecting the first dot to 
the lead brings this value to G = 1.98 Go- 



V. CONCLUSION 

In conclusion, we have demonstrated how to tune a 
linear array of quantum dots coupled to superconductors 
in presence of Zeeman field and spin-orbit coupling to 
resemble the Kitaev chain that hosts Majorana bound 
states at its ends. Furthermore, we have presented a de- 
tailed procedure by which the system is brought to the 
so-called "sweet spot" in parameter space, where the Ma- 
jorana bound states are the most localized. This proce- 
dure involves varying the gates potentials and supercon- 
ducting phases, as well as monitoring of the excitation 
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FIG. 4: Conductance as a function of bias voltage during 
simultaneous tuning of G± and V n for the first pair of dots. 
The three different plots represent the situation before (dotted 
line), at an intermediate stage (dashed line), and after (solid 
line) the tuning. The arrow indicates the evolution of the first 
peak upon tuning, and the number of successive changes of 
G± and V n are shown for each curve. By bringing the first 
peak to zero, the third tuning step is achieved. 
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FIG. 5: Conductance as a function of bias voltage for a sys- 
tem composed of three tuned quantum dots (dashed line). 
The zero bias peak signals the presence of Majorana bound 
states at the ends of the chain. The first and second excited 
states are consistent with those expected for a three-site Ki- 
taev chain, namely E\ = 2t\ and E2 = 2fe (vertical dashed 
lines), given the measured values of t\ — Ai and £2 = A2, 
obtained after finalizing the two dot tuning process. As de- 
scribed in the main text, after increasing the transparency 
of the lead QPC leads to a zero bias peak having a height 
G = 1.98G (solid line). 



spectrum of the system by means of resonant Andreev 
conductance. 

We have tested our procedure using numerical simula- 
tions of a system of three QDs, defined in a 2DEG, and 
found that it works in systems with experimentally reach- 
able parameters. It can be also applied to systems where 
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FIG. 6: Geometry of the quantum dot chain. The quantum 
dots have a width Wdot and length equal to Ldot- Quantum 
point contacts have a longitudinal size Lqpc and a transverse 
dimension equal to either Ldot or Wuot- Leads are semi- 
infinite in the x direction, and superconductors are modeled 
as semi-infinite systems in the y direction. 



quantum dots are denned by other means, for example 
formed in a one-dimensional InAs or InSb wire. 
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Appendix A: System parameters in numerical 
simulations 



In this section, we describe the parameters used 
throughout the numerical simulations. The quantum 
dots and quantum point contacts are modeled using a 
tight-binding model defined on a square lattice, with 
leads and superconductors taken as semi-infinite. 

The characteristic length and energy scales of this sys- 
tem are the spin-orbit length /go — h 2 /ma, and the spin- 
orbit energy Ego — ma 2 /h 2 . We simulate an InAs sys- 
tem in which the effective electron mass is to = 0.015m e , 
where m e is the bare electron mass, taking values of 
E so = 1 K = 86 [ieV and l so = 250 nm. 

We consider a setup composed of three quantum dots, 
like the one shown in Fig. [6} Each of the three dots has a 
length of £dot = 208 nm and a width PFdot = 104 nm. 



Quantum point contacts have a longitudinal dimension 
of Lqpc = 42 nm, which is the same as the Fermi wave- 
length at quarter filling. 

The value of the hopping integral becomes t = 
h 2 /(2ma 2 ) — 55.8 meV, with a = 7nm. Disorder 
is introduced in the form of random uncorrelated on- 
site potential fluctuations, leading to a mean free path 
^mfp = 218.8 nm. The system is placed in a perpendic- 
ular magnetic field characterized by a Zeeman splitting 
V z = 336^eV, which, given a g-factor of 35K/T, cor- 
responds to a magnetic field B z — lllmT. Each dot is 
additionally connected to a superconductor characterized 
by a pairing potential |A SC | = 0.86 meV. 

The potential profile across a quantum point contact 
is given by 



Vqpc(;e) 



tanh 



x - 



tanh 



w 

L V 2 

w 
~ 2 



(Al) 



where x 6 [—L/2,L/2] is the transverse coordinate across 
the quantum point contact, h is the maximum height of 
^QPd s" fixes the slope at which the potential changes, 
and w is used to tune the QPC transparency. Two ex- 
amples of potential profiles are shown in Fig. [7] 




x [nm] 



FIG. 7: Potential profile Vqpc(:e) across the transverse direc- 
tion of a quantum point contact. For the maximum value of 
this potential, no states are available for quasiparticles in the 
2DEG. The two curves show potential profiles for two different 
QPC transparencies, corresponding to s = 17 and w = 87.4, 
39.5 nm for the solid and dashed curves respectively. 
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